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Abstract. Several methods for the determination of the mass of the top quark with 
the ATLAS detector at the LHC are presented. All dominant decay channels of the 
top quark can be explored. The measurements are in most cases dominated by system- 
atic uncertainties. New methods have been developed to control those related to the 
detector. The results indicate that a total error on the top mass at the level of 1 GeV 
should be achievable. 
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1 Introduction 

A precise measurement of the top quark mass will be a main goal of top physics 
at the LHC. The combined top mass value from the Tevatron run I is m,t = 
174.3 ±5.1 CeV H], and the expected accuracy obtained after run H will be 3 
GeV 0. 
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Motivations for an accurate determination of the top quark mass are nu- 
merous. It is a fundamental parameter of the Standard Model (SM) and should 
therefore be measured as precisely as possible. An accurate value of the top mass 
would help to provide a rigorous consistency check of the SM and to constrain 
some parameters of the model such as the mass of the Higgs boson. Further- 
more, a high level of accuracy on the top mass value (for example improving 
the accuracy down to Arrit ^ 1 GeV) is also desirable, both within the SM and 
the Minimal Supersymetric Standard Model (MSSM) framework T. In the SM, 
such an accuracy would significantly improve the precision on the W boson mass 
prediction while in the MSSM, it would put constraints on the parameters of 
the scalar top sector and would therefore allow sensitive test of the model by 
comparing predictions with direct observations. 

Because the top quark, as other quarks, cannot be observed as a free particle, 
the top quark mass is a purely theoretical notion and depends on the concept 
adopted for its definition. With increasingly-precise measurements on the hori- 
zon, it is important to have a firm grasp of exactly what is meant by the top 
quark mass. Thus far the top quark mass has been experimentally defined by 
the position of the peak in the invariant mass distribution of the top quark's 
decay products, a W boson and a b quark jet. This closely corresponds to the 
pole mass of the top quark, defined as the real part of the pole in the top quark 
propagator. 

The renormalisation scale dependence is less than 10 MeV for the range of 
the scale between 30-150 GeV. The dominant theoretical uncertainty for the top 
mass caused by uncertainty on the strong coupling constant is less than 150 
MeV to the binding energy, which would give uncertainty of 75 MeV in the pole 
mass. Corresponding theoretical uncertainty in the MS mass would be about 
±12 MeV(4j. However, this definition is still adequate for the analysis of top 
quark production at LHC where uncertainty in the top mass measurement will 
be of order 1 GeV. Because of fragmentation effects it is believed that the top 
quark mass determination in an hadronic environment is inherently uncertain to 
0{Aqcv) 00. 

At the LHC, the top quark will be produced mainly in pairs through the 
hard process gg tt (90% of the total tt cross-section) and qq ti (remaining 
10% of the cross-section). The next-to- leading order cross-section prediction for 
tt production is a{tt) = 833 pb 0. Thus the LHC will be a top factory as more 
than 8 million tt pairs will be produced per year at low luminosity (corresponding 
to an integrated luminosity of 10 fb~^). The electroweak single top production 
processes, whose cross-sections are in total approximately one third of those of 
tt production, have not been investigated for the determination of the top mass. 

Within the SM, the top quark decays almost exclusively into a W boson and 
a b-quark {t Wb). Depending on the decay mode of the W bosons the events 
can be classified into three channels: the lepton plus jets channel, the dilepton 
channel and the all jets channel. In the lepton plus jets channel, one of the W 
boson decays leptonically {W Iv) and the other one hadronically {W jj). 
Considering electrons and muons, the branching ratio is BR = 2x2/9x6/9~ 
30%. The final state topology is gg tt ^ {ijb){lvh). In the dilepton channel, 
both W bosons decay leptonically with BR = 2/9 x 2/9 ~ 5%. In the all jets 
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channel, both W bosons decay hadronically with BR = 6/9 x 6/9 ~ 44%. 

This paper summarizes studies of the top mass measurement, including up- 
dates of studies presented in the ATLAS Technical Design Report ,7| as well as 
several new analysis. 

Unless otherwise indicated, all analysis were performed with events generated 
with Pythia 8 and passed through the ATLAS fast simulation package Atlfast 
|H] for particles and jets reconstruction and momenta smearing. Jets are defined 
as massless objects by summing the momenta of clusters of energy deposited in 
the calorimeters. Clusters are associated to form jets using a cone algorithm with 
S{R) < 0.4. A tagging efficiency of 60% for b-jets was assumed. Cross-checks of 
some results have been made using the detailed GEANT-based simulation of the 
ATLAS detector. The top mass is extracted by an analytic fit to the event by 
event reconstructed invariant mass. 

2 Top mass measurement in the lepton plus jets channel 

The lepton plus jets channel is probably the most promising channel for an 
accurate determination of the top quark mass. Three methods to measure the 
top mass are envisaged. The simplest method consists in extracting the top 
mass from the three jets invariant mass of the hadronic top decay. In the second 
method, the entire tt system is fully exploited to determine the top quark mass 
from a kinematic fit. In the last method, the top mass is still determined from a 
kinematic fit, but the jets are reconstructed using a continuous algorithm. 

2.1 Event selection and background rejection 

Taking into account the total tt cross-section and the branching ratio, one can 
expect 2.5 millions ti pairs with this topology to be produced per year assuming 
an integrated luminosity 10/5^^. 

The signal final state ti — ^ Wb Wb jjb Ivb (with I — e, n) is characterized 
by one high transverse momentum lepton, large transverse missing energy iJ^""^, 
and high jet multiplicity. The following background processes were considered: 
bb lv + jets, W + jets — > lv + jets, Z + jets +jets, WW Iv + jets, 

WZ ^ Iv + jets, and ZZ l^l^ + jets. At production level, the signal over 
background ratio is very unfavorable (S/B ^ 10^'''). 

A high level of rejection was obtained using the following requirements: one 
isolated lepton with px > 20 GeV, i;™*^** > 20 GeV, and at least four jets 
reconstructed with a cone size of AR — 0.4 with pT > 40 GeV and \t]\ < 2.5, 
of which at least two are tagged as b-jets. The efficiency of the selection for the 
various background processes is shown in tabled After selection cuts, the signal 
over background ratio is extremely good (S/B ^ 78), and the remaining number 
of signal events is approximately 87000 (for an integrated luminosity of 10~^fb). 

The requirement of having at least two b-tagged jets in the final state helps in 
rejecting a large part of the physical background, but also reduces considerably 
the signal sample. The fraction of signal events with at least two b-tagged jets 
is three times smaller than the fraction with at least one b-tagged jet (see figure 

Requiring only one b-tagged jet would decrease the signal over background 
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Process 


Cross-section 


Total efficiency 




(Pb) 


(%) 


tt signal 


250 


3.5 


bb liy + jets 


2.2 X 10^ 


3 X 10-*^ 


W + jets ^11^ + jets 


7.8 X 10^ 


2 X 10--* 


Z + jets ^ l+r + jets 


1.2 X 10^ 


6 X 10-5 


WW ^lv + jets 


17.1 


7 X 10-3 


WZ -^lu + jets 


3.4 


1 X 10-2 


ZZ l+l- + jets 


9.2 


3 X 10-3 



Table 1. Cross-section and selection efficiency for signal and background pro- 
cesses. For the background, the hard scattering processes are generated with a 
cut on the transverse momentum at 20 GeV 
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Fig. 1. tt sample as a function of the number of b-tagged jets in the event. 



ratio from 78 to 28, which would still be acceptable. This extended sample can 
be used for the top mass measurement using the hadronic top decay, requiring 
the tagged b-jet to be the one belonging to the hadronic top final state. 

2.2 Top mass measurement using the hadronic top decay 

In the following method, the top mass will be determined from the invariant 
mass of the three jets arising from the hadronic top decay {t — * Wb jjb). To 
accurately reconstruct the decay, one should: i) identify the jets associated to the 
hadronic top decay among all other jets, ii) precisely calibrate the jet energies 
and directions. 

2.2.1 Jet association At least four jets are expected in the event: two from 
the hadronic W decay and two b-jets. Additional jets will be produced by initial 
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Fig. 2. Dijet invariant mass distributions for events with at least 1 b-tagged jet 
(left plot), and at least 2 b-tagged jets (right plot). The shaded area represents 
the combinatorial background. Both plots are for 10 fb~^. 



state radiation (ISR) and final state radiation (FSR) effects. The association of 
jets to the original partons is done as follows. The hadronic W decay is first 
reconstructed: all the non b-tagged jets are paired together and the jet pair with 
an invariant mass closest to the W mass is taken as the right combination for 
the W. The di-jet invariant mass distributions for events selected by requiring 
two b-tagged jets or at least one are shown on figure|2 When the two associated 
jets are reconstructed, 80 % of the true W decays are selected, which is realized 
in 45 % of the cases. This leads, in a mass window \Mjj - Af^^'^l < 20 GeV, 
to a purity of 55 % for events selected with at least one b-tagged jet and 66 % 
for events selected with two b-tagged jets. The width of the mass distribution is 
7.4 GeV for both cases. 

The next step is to associate a b-tagged jet to the reconstructed W. When 
events are selected with only one b-tagged jet, the association is performed if 
the b-tagged jet is closer to the reconstructed W than to the isolated lepton 
{AR{b,W) < AR{b, lepton)). The efficiency of this association criteria is 82%. 
In the presence of two b-tagged jets, the chosen b-tagged jet is the one giving 
the highest pT for the reconstructed top, giving an association efficiency of 81%. 
The reconstructed three jets mass distributions are presented on figureOlfor the 
two event selection criteria. The peak width is 12 GeV in both cases. 

The overall association purity and efficiency within a top mass window of 
±35 GeV around the top mass peak are summarized in table |21 The top mass 
determination will not be limited by the statistics even if the analysis is restricted 
to the two b-tagged jets sample. Nevertheless, the large one b-tagged jet sample 
allows further dedicated cuts with negligible impact on the statistical precision 
of the top mass determination. In the sample with two b-tagged jets, the overall 
reconstruction efficiency is 1.2%, leading to 30000 events for one year of running 
at low luminosity (per 10 fb~^). For simplicity, only this case will considered in 



Investigation of top mass measurements with the ATLAS detector at LHC 



7 




100 200 300 400 



M..(GeV) 



> 




100 200 300 400 



M,.,(GeV) 



Fig. 3. Top mass distributions for the 1 b-tagged jet sample (left plot) and the 
2 b-tagged jets sample (right plot). The shaded area represents the combinatorial 
background. Both plots are for 10 fb~^ . 



the following. 

1 b-tagged jet sample 2 b-tagged jets sample 
Top purity (%) 65 69 
Total efficiency (%) 2^5 1^2 

Table 2. Summary table for the two samples considered. Purity and total effi- 
ciency are related to a top mass window of ±35 GeV around the generated mt 
value. 



2.2.2 In situ jet energy and direction calibration As the top mass is 
determined from the invariant mass of a three-jet system, the accuracy of the 
measurement depends on how well the jet energies and directions are recon- 
structed. A mis-measurement of 1% of the jet energies induces a top mass shift 
of 1.6 GeV. Similarly, a mis-measurement of 1% of the cosinus of opening angle 
between the two jets from the W and between the reconstructed W and the b-jet 
induces a top mass shift of 1.2 GeV. Therefore an excellent absolute energy scale 
and angle measurement are required to precisely determine the top quark mass. 

Numerous effects have to be taken into account to determine the initial par- 
ton energy from the energies deposited in the ATLAS calorimeters. Prior to 
data taking, the accumulated knowledge on the detector performances and char- 
acteristics, on physics effects like initial and final state radiations, underlying or 
minimum bias events plus the impact of the jet finding algorithms will allow to 
reach a 5-10 % on the absolute energy scale '7'. In-situ calibrations will fix the 
absolute energy scale through the study of known processes, taking into account 
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in a global way the remaining inaccuracies on the knowledge on the various 
effects described previously. 

It has been shown that an accurate absolute energy calibration of light quark 
jets and b-jets can be extracted from Z+jct events [71 llUj. within an expected 
precision of about 1 %. However, this calibration applied to the W mass recon- 
struction 11, 12 , leads to a shifted W mass. Due to the energy sharing between 
jets, the opening angles are systematically underestimated leading to a less pre- 
cise mass measurement 

Here, to avoid any dependency from external inputs, it is proposed to perform 
an situ calibration in which both the absolute energy and direction calibration 
are extracted from the W jj channel itself. For this purpose, a cleaner sample 
of W candidates has been selected from the tt events. Initially, the jets are 
not calibrated but corrected for cell energy sharing effects. In addition to the 
preselection cuts, the di-jet invariant mass is required to fall within a mass 
window of ± 20 GeV around the peak value and the three-jet invariant mass to 
fall within a mass window of ± 15 GeV around the peak value. 

The in-situ calibration is performed through a minimization procedure in 
which the dijet mass is constrained to the known W mass. 

Non-calibrated jet energies are shifted from the initial parton energies due 
mainly to the jet cone algorithm procedure and FSR effects. However the induced 
shift is in general smaller than the energy resolution of the jets. This allows to fix 
the term to the intrinsic calorimeter energy resolution. The same approach 
is employed for the jet directions. 




Fig. 4. Calibration factors obtained 

, , , J J. ■ J Fie. 5. Final jib invariant mass distri- 

event by event, and parametrization aj- ' 

, 7 1, J- n, bution. 
ter calibration fit. 

An energy correction factor K is obtained at the end of the fitting procedure, 
for each jet and for each event. The distribution as a function of the raw initial 
jet energy is shown in figure ^ Finally, the function K = P2 + ^raw is fitted 
to the distribution (the resulting parameters are shown on the figure) leading to 
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Fig. 6. Results of the calibration fit on jet energy and direction. The two solid 
lines show the ±1% precision level. 



a calibration function, without an a priori knowledge on the initial calibration 
function shape. In figure El the comparison between the initial parton energy 
and the calibrated jet energy at various steps of the procedure is represented. 
One can see the impressive effect of the in-situ calibration procedure as the ratio 
of the parton energy to the reconstructed jet energy remains well below the 
level of 1%. The improvement brought by this procedure is also clear on the 
reconstruction of the opening angle between the two jets, as seen on figure El It 
should be noted that the combinatorial background does not introduce a sizable 
bias on the calibration factors. 



2.2.3 Top mass reconstruction The selected three jets invariant mass dis- 
tribution is shown in figure El The light quark jets were calibrated as described 
above and the b-quark jets were calibrated using Z+b events [Jj. The mass peak 
is in agreement within 100 MeV of the generated value. The peak width is around 
11 GeV, leading to a statistical error on the top mass of the order of 100 MeV 
for one year of running at low luminosity (per 10 fb~^). 

2.2.4 Full simulation results The analysis presented in section 2.2 of the 
determination of the top quark mass using the hadronic top decay has been 
repeated using fully simulated events. For this purpose, 30000 tt events were 
processed through the GEANT-based ATLAS detector simulation package. 

Events were generated under restrictive conditions at generator level. These 
conditions include, for example, cuts on the transverse momentum of the tt decay 
products. Therefore any direct comparison with the results presented in section 
2.2 should be avoided. The comparison is made using the same generated events 
which have been passed through both the fast and full simulation package. 

FigurejZlrepresents the mjj and rrijjb distributions for fast and full simulation. 
The rrijj invariant mass resolution is 7.3 GeV for fast simulation and 8.1 GeV for 
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Fig. 7. Top mass determination using the hadronic top decay: rrijj and rrijjb 
invariant mass distributions. The left-handed plots represent the distributions 
obtained from fast simulation and the right-handed plots represent the distribu- 
tion obtained from full simulation. 
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full simulation. The rrijji, invariant mass resolution is 11.4 GeV for fast simulation 
and 13.4 GeV for full simulation. In the top mass window 175 ± 35 GeV, the 
signal purity and overall efficiency are P = (79 ± 2)% and E = (6.4 ± 0.5)% for 
fast simulation and P (78 ± 2)% and E = (5.7 ± 0.5)% for full simulation. 
These results are summarized in table |3| 



Quantity 


Fast simulation 


Full simulation 


nijj resolution (GeV) 


7.3 


8.1 


rujjb resolution (GeV) 


11.4 


13.4 


Signal purity (%) in 175 ± 35 GeV 


79 


78 


Signal efficiency (%) in 175 ± 35 GeV 


6.4 


5.7 



Table 3. Top mass determination using the hadronic top decay: comparison 
between fast and full simulation. 



The results obtained for the signal purity and overall efficiency as well for the 
rrijj and mjjb invariant mass resolutions are in reasonable agreement between 
fast and full simulation, though the resolutions from GEANT are somewhat 
worse. In addition, the shape and amount of the combinatorial background for 
both the W and top masses reconstruction are also in good agreement between 
the two types of simulations. 

2.2.5 Systematic uncertainties To estimate the effect of an absolute jet 
energy scale uncertainty, different miscalibration coefficients were applied to the 
reconstructed jet energy. A top mass shift per percent of miscalibration was 
obtained. For light quark jets, the effect is small as the jet are re-calibrated 
in-situ. For b-quark jets, a 1% miscalibration induces a top mass shift of 0.7 
GeV. 

The presence of initial state radiation of incoming partons (ISR) and final 
state radiation from the top decay products (FSR) can impact the measure- 
ment of the top mass. To estimate their effect, a top mass shift due to ISR was 
computed as the difference between the value of the top mass determined with 
ISR switched on (usual data set) and ISR switched off. The same approach was 
employed for FSR. The level of knowledge of ISR and FSR is of order of 10%. 
Therefore, as more conservative estimate, the systematic uncertainty on the top 
mass was taken to be 20% of the corresponding mass shifts. 

The b-quark fragmentation was described by the Peterson fragmentation 
function [14, . This function is parametrized in terms of one variable £(,. The 
default value was set at et — —0.006, with an uncertainty of 0.0025 ^S]. The 
top mass was determined with another sample of events generated with Sb = 
—0.006 -I- 0.0025. The difference of the top mass value between the usual sample 
and the latter was taken to be the systematic uncertainty on the top mass due 
to the knowledge of £{,. 

Uncertainties due to the combinatorial background (which is the main back- 
ground) were also estimated by varying the assumptions of the background shape 
and size in the fitting procedure. Fits of the three-jet invariant mass distribution 
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Systematics 



Suit (GeV) 



Light jet energy scale 
b jet energy scale 



0.2 

0.7xx% 



Initial State Radiation 
Final State Radiation 
b-quark fragmentation 



0.1 
1 

0.1 

0.1 
0.1 



Combinatorial background 
Statistical error 



Table 4. Summary of systematic errors in the inclusive lepton plus jets sample. 
X represents the level at which the b jet energy scale will be known. 



were performed using a Gaussian shape for the signal and either a polynomial 
or a threshold function for the background. The resulting systematic error on 
the top mass was 0.1 GeV. 

All the results are summarized in tabled The main contributions are from 
FSR and b-quark energy scale. Adding in quadrature all the contribution leads to 
a total systematic uncertainty on the top mass of the order of 1.3 GeV, provided 
the b-quark jets can be calibrated within 1%. 

2.3 Top mass measurement using a kinematic fit 

In the previous section, it was shown that the top quark mass can be mea- 
sured in the lepton plus jets channel with an accuracy better than 2 GeV. This 
error is totally dominated by systematic effects, in particular the b-quark jet 
energy scale and FSR. In order to reduce further the systematic uncertainties, 
another method is proposed in the following, where the entire tt final state is 
reconstructed by a kinematic fit. This method is aimed to reduce the impact 
of poorly reconstructed jets (due to effects arising from FSR and semi-leptonic 
decays of b-quark jets). The final state can divided into two parts: i) the leptonic 
part corresponding to the leptonic top decay {t{t) — > lvb{b)) and ii) the hadronic 
part corresponding to the hadronic top decay {t{t) jjb{b)). 

The hadronic part is reconstructed in a similar way to the previous section. 
For the light quark jets, the absolute energy scale was taken from the in-situ 
calibration described above and for b-quark jets, it was taken from Z+h events 
|7]. The invariant mass distribution of the selected di-jet pairs is reconstructed 
with a width on 6.2 GeV (see figure O. The three-jet system is reconstructed 
with a width of 11.2 GeV (see figure |SJ). Over the entire mass range, the purity 
is 51% and is increased to 71% within a mass window of ± 35 GeV around the 
peak value. 

The leptonic final state cannot be directly reconstructed due to the presence 
of the undetected neutrino. Nevertheless, the neutrino four-momentum can be 
estimated in two steps. First, the transverse component of the neutrino mo- 
mentum can be approximated by the transverse missing energy (see figure . 
The longitudinal component of the neutrino momentum can then be deduced 
with a quadratic ambiguity, by constraining the invariant mass of the lepton- 
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Fig. 8. Hadronic part: reconstructed W mass (left plot) and reconstructed top 
mass (right plot). The combinatorial background contributions are shown (shaded 
areas) together with the full invariant mass distributions (full line). 



> 




Fig. 9. Ratio of the neutrino trams- 
verse momentum at parton level over 
the reconstructed neutrino transverse 
momentum. 



Fig. 10. Reconstructed top mass from 
the leptonic part. The combinato- 
rial background contribution is shown 
(shaded area) together with the full in- 
variant mass distribution (full line). 



neutrino system to the known W mass value. Finally, the remaining b-tagged 
jet is associated to the reconstructed W. In most of the cases, there are only two 
b-tagged jets present in the event, and one has already been associated. In case 
of additional b-jets, the closest one to the isolated lepton is chosen. Per event, 
two leptonic top masses are computed, corresponding to the two neutrino pz 
solutions. The distribution of the leptonic top mass the closest to the hadronic 
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Fig. 11. Le^t plot: distribution. Right plot: top mass after kinematic fit. 



mass is represented in figure ^1 It is fitted by a third order polynomial plus a 
Gaussian, leading to a peak width of 12.4 GeV, similar to the hadronic resolu- 
tion. The event is eventually kept if one of the two leptonic top masses is within 
the same top mass window as the hadronic part. 

Therefore, the entire ti final state can be reconstructed, with a twofold am- 
biguity due to the neutrino reconstruction. After selection and mass window 
cuts, the final sample is composed of 18000 signal events and 7000 combinato- 
rial events (the contribution from other physical background processes is totally 
negligible). The total signal {tt events with all jets well assigned) efficiency is 
0.7% with a purity of 73%. 

2.3.1 Top mass determination The kinematic fit is performed in such a 
way that the jets and lepton energy, the jets direction (in terms of rj and </>) and 
the three components of the reconstructed neutrino momentum can vary freely 
within their corresponding resolutions. The following kinematic constraints were 
employed: 

mjj = M^^^, mi^ = M^DG ^-^..^ ^ ^^^^^ ^ ^fit^ 
On an event-by-event basis and for both neutrino solutions, a is minimized 
[T^ . The output of the fit is the top mass estimator Mf^p. The neutrino solution 
with the lowest is selected. 

The shapes of the distributions are different for signal and combinatorial 
background, as is shown in figure ^2 a cut at < 4 increases the tt purity to 
more than 80%. The fitted top mass is also shown in figure ITTI 

The constraints which define the are strong on the di-jet and the lepton- 
neutrino systems, due to the good knowledge of the W mass, but give a poor 
relative constraint on the two b-tagged jets. As a consequence, the x^ is directly 
related to the quality of the reconstruction of the b-jets, particularly as the b-jet 
energy can be underestimated when FSR or leptonic decays occurs. Furthermore, 
the quality of the top mass reconstruction can be altered, with a mass value un- 
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Fig. 12. Left plot: relative difference between b-jets and b-partons energies. 
The shaded area defines the "badly reconstructed" b-jets. Right plot: probability 
to have the two b-tagged jets "well reconstructed" as a function of ■ 



derestimated in correlation with the effects on the b-jets energy reconstruction. 
The left plot of figure El shows the relative difference between the b-jets and 
b-partons energies. The jets which belong to the tail of this distribution can be 
defined as badly reconstructed jets. The probability to have the two b-jets "well 
reconstructed" decreases when the increases (see figure [T^ . The top mass 
dependence on the is shown in figure El It can be noticed that the fitted top 
mass becomes independent on the x^ if only events with "well reconstructed" 
b-jets are kept. To some extend, the x^ value allows to distinguish between "well 
reconstructed" b-jets from the others. 

The top mass is estimated in the following way. Equal samples per slices 
of are built. The top mass is computed for each sample, from a Gaussian fit 
around zbl.Scr of the mass peak. Figure IT!^ presents two examples, for slices of 
with mean values x^ — 0.12 and x^ = 3.63. Finally, the top mass is determined 
as mt = MjQp(x^ = 0), from a fit by a linear function to the distribution (see 
figure EJ- III one year of running at low luminosity (per 10/5^^) the statistical 
error would be 120 MeV. 

2.3.2 Systematic uncertainties The study of the systematic uncertainties 
were handled in the same way as in the previous section. The results are pre- 
sented in table El For FSR, taking 20% of the mass shift obtained between FSR 
switched on and off leads to a systematic error on the top mass of 0.5 GeV. How- 
ever, this estimate is an upper limit as the mass shift takes into account effects 
due to a wrong b-quark jet calibration. This effect could possibly be reduced: 
if the absolute b-jet scale is obtained with Z+h events generated with FSR off, 
the systematic error would be decreased to 0.1 GeV. This is not surprising, since 
when the size of the FSR contribution increases, only the slope in the figurelTHis 
modified. This demonstrate that events with large FSR contributions populate 
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Fig. 13. Top mass distributions for different values of . One can observe that 
for the highest the corresponding top mass is lower, and the mass resolution 
is higher. The combinatorial background contributions is shown (shaded area) 
together with the full invariant mass distribution (full line). 
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Fig. 15. Reconstructed top mass ver- 
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high values. 

The Hnearity of the method has been checked using several samples generated 
with different top masses. The same algorithm was used on all samples, with the 
same calibration functions for light and b jets. As is shown in figure El the 
estimated top mass depends linearly on the generated top mass. 

In table the errors are divided into three types: i) the statistical error, ii) 
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Internal Systcmatics 


Suit (GeV) 


light jet energy scale 


0.2 


Initial State Radiation 


0.1 


Final State Radiation 


< 0.5 


b-quark fragmentation 


0.1 


Combinatorial background 


0.1 


Total 


< 0.6 


Statistical error 


0.1 


b jet energy scale 


0.7xa:% 



Table 5. Summary of the systematic errors for the top mass in the inclusive 
lepton plus jets sample, when the top mass is reconstructed using a kinematic fit. 



internal systematic errors and iii) the systematic error due to the b-quark jet 
absolute energy scale which depends on external inputs. No systematic errors 
have been accounted for the Monte Carlo description of the decay as it has been 
demonstrated to be negligible (I2J. Assuming the b-quark jet absolute energy 
scale can be determined within 1%, the total error on the top mass measurement 
is of the order of 0.9 GeV. 
In summary: 

mt = m ± 0.1 (stat.) ± (0.3 - 0.6) (internal syst.) (±0-7 X x) (external syst.) 

where x — b-quark jet miscalibration in %. 

2.4 An additional technique: mass measurement using a continuous 
jet algorithm 

The main sources of systematic uncertainty entering in the top mass measure- 
ment in the inclusive lepton plus jets channel arise from the non-precise knowl- 
edge of the correction factors applied to the jet energy and of some physics 
processes parameters. In order to reduce the impact of the latter effects (mainly 
final state radiation) in the top mass determination, an approach based on the 
continuous definition of jets has been investigated |16| . 

This approach is based on the following ingredient: a continuous jet definition 
and mass estimation. This idea was first introduced in ^^l- It is based on the 
consideration that the discrete nature of the standard jet definition may cause 
some problems. Two of these problems are: i) from the mathematical point of 
view, the breakdown of continuous distributions gives rise to instabilities (large 
statistical fluctuations) and therefore increases the statistical error on the mea- 
surement result, ii) the transition from a continuous energy deposition in the 
calorimeter to a fixed structure jet causes the loss of important information. As 
an example, let us consider the case of final state radiation. 

As a result, the fraction of the total energy contained in the cone around 
the quark direction in any standard jet finding algorithm is subject to large 
fluctuations which cannot be precisely predicted from theory. This induces a large 
contribution to the systematic errors in the mass measurement. A continuous jet 
definition allows instead to reduce the dependence of the estimation of the details 
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of the jet shape description in the Monte Carlo and in the jet reconstruction 
procedure. 

In addition, this method can be coupled to a jet energy calibration based on 
the W mass measurement and on the reconstruction of the tt final state with a 
constraint kinematic fit, as presented in the previous section. 

2.4.1 Method The event selection criteria are similar to the ones used in the 
previously described analysis (see section 2.1). The continuous jet algorithm has 
been realized in the following way. Initially, a "fixed cone" jet finding algorithm is 
used, with a definite cone size. For the same event, the analysis is then repeated 
varying the jet cone size. Here, the cone sizes range from 0.3 to 1.0 with step of 
0.1. 

In a preliminary stage of the analysis for each cone size a jet energy correction 
factor has been defined by calculating the two-jet invariant mass distribution for 
non b-tagged jets with pr > 40 GeV. The W-peak position has been fitted using 
the sum of a Gaussian and Tchebyshev polynomials up to the fourth order. 

The same corrections as for the light quark jets have been applied for the 
b-jets, even if not being the optimal ones. Two non b-tagged jets with px > 40 
GeV have been selected for which the combined invariant mass was close to 
the mass of the W peak for a given cone size chosen for the jet reconstruction 
algorithm. The requirement |M — Mwtr-^^ \ < 25 GeV has been used. A b-quark 
jet energy rescaling has been applied to the b-tagged jets with pt > 40 GeV, 
according to the correction factor for the given cone size. At least two b-quarks 
are required in each event. If there are more, the combination with the highest 
Pt is selected. 

For the given cone size, a constraint fit procedure was then applied to the 
selected jets. This kinematic fit uses three constraints rrijj — Mw, mii, = Mw 
and mjjyf, = "^jjb which, together with the missing momentum measurement, 
allows to determine the neutrino momentum unambiguously. Only the energies 
of the jets are tuned during the fit. However, in the method proposed here, a 
slightly different approach has been exploited to improve the accuracy of the 
reconstruction procedure: a robust modification |18lll9j of the fitting functional 
have been used instead of the usual J^l ■ 

Only the events with a small minimum value of the fitting functional have 
been retained. This value is not a because of the robust fitting method and 
depends on jet and missing momentum error parameterizations. The invariant 
mass of the 3-jet system obtained after the constraint fit has then been considered 
as the top quark mass estimate for the current event and current jet cone size. 

After having a top quark mass estimate for a certain jet cone size, the cone 
size is changed and the whole procedure is repeated for the given event, starting 
again from the jets finding. If the jet finding procedure is unable to find the 
required minimum number of jets with pT > 40 GeV the algorithm skips to the 
next event. All the top quark masses calculated for all the events and all the 
jet cone sizes are finally summed up in one single histogram. The distribution 
obtained as the result of the reconstruction procedure is shown in the left plot 
of figure El 

For different events, the mass value at the peak position is obtained for 
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Fig. 16. Left plot: Summary distributions of the top quark mass estimates for 
all cone sizes. Right plot: Difference between the generated and reconstructed 
top quark masses, for different top mass values, after having applied the mass 
reconstruction procedure as described in the text. The result of a linear (constant) 
fit is also shown. 



difFerent cone sizes. The point here is that a position of accumulation should be 
a more robust estimator of the jet system invariant mass with respect to the 
single value which is obtained by calculation of only one mass for each event. In 
figure rm the 3-jet system invariant masses (the top quark mass estimator) are 
shown as a function of the cone size AR as well. The summed distribution for 
all the AR values is shown at the bottom of figure [T7I 

The top quark decay peak is then fitted with a Breit-Wigner function includ- 
ing a 4th order Tchebyshev polynomial to describe the background. The position 
of the peak is considered as the final estimate of the top quark mass. The Breit- 
Wigner shape has been selected because it gives a much better description of 
the signal compared with a Gaussian shape. However, an equally good may 
be obtained by using a sum of two Gaussian with the same mean for the signal 
description. The difference between the results with the two methods has been 
treated as a systematic error. 

2.4.2 Results A typical mass distribution obtained by applying the top quark 
mass reconstruction procedure described above is shown in figure IT^ for 400000 
tt generated events with mf^'-^ = 175 GeV. 

A mean statistical error of the top quark mass estimation obtained with the 
fit of the mass distribution in figure [THl bv a Gaussian function plus a 4th order 
polynomial background is Smt ^ 100 MeV, for 400000 generated events. 

In order to evaluate the statistical properties of the top quark mass estima- 
tion procedure from Smt^ some additional steps are needed. Due to the contin- 
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Fig. 17. Mass distributions for the three jets system and for various cone sizes. 
The vertical line indicates the generated top mass value. The bottom-right plot 
represents the sum of all distributions. 



uous jet definition method and to the data based b-jet energy correction, the 
statistical error given by the fit of the distribution is not correct. The invari- 
ant mass distributions obtained from the same event sample with different cone 
sizes used by the algorithm are strongly correlated and cannot be treated as 
independent. The statistical error of the b-jet energy correction factors has to 
be taken into accoimt. Both effects lead to an imderestimate of the statistical 
error as obtained from the invariant mass distribution fitting procedure. For the 
correction, one should rescale statistical errors from the fitting procedure. On 
the other side the mass distributions obtained with the different jet cone sizes 
are not the same because even the number of reconstructed jets in event may 
be different for AR = 0.3 and AR = 1.0, then the summed distribution has 
more information than one obtained with the single AR. So one must rescale a 
statistical error obtained in the fitting procedure by a factor in the range l...-\/8. 
The determination of such rescaling factor is not an easy task to be performed 
analytically. However, one can easily determine the corresponding correction in 
Monte Carlo experiment, with the help of pull distribution. For each of the top 
quark masses {mt = 165, 170, 175, 180, 185 GeV), 400000 events were generated 
and reconstructed. This procedure has been repeated five times. The differences 
between the generated and reconstructed top masses are determined for all input 
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top masses and the corresponding pull distribution is obtained. The dispersion 
of the pull distribution (equal to the scaling factor) is 1.6. After multiplying 5mt 
by 1.6, taking into account b-tagging efficiency and rescaling to an integrated 
luminosity of 10 fb^^, an estimated statistical error about Smt ~ 65 MeV is 
found. The final statistical error of the top mass will be less than 100 MeV. 

2.4.3 Systematic errors As expected, and as can be seen in figure El the 
reconstructed three-jet mass depends on the jet cone size. The final value ex- 
tracted for the top mass will, therefore, depend on the choices for the minimum 
and maximum cone sizes to use in the analysis. The analysis, which used a cone 
size step of 0.1 with minimum and maximum values of 0.3 and 1.0 respectively, 
was re-done excluding either the minimum or maximum values and again by 
shifting the grid of cone sizes by 0.05. The largest change observed in the resul- 
tant top mass, namely 250 MeV, has been assigned as the systematic error in 
the top mass due to uncertainty in the range of cone sizes to use. 

The effect of the cut value on the determination of the top mass has been 
checked. The difference between the reconstructed top mass and the generated 
value has been plotted versus the cut value. A maximum difference of 200 
MeV was found. This value has been taken as the systematic error due to the 
X^ dependence of the top quark mass determination. 

The description of the background and signal shapes in the fitted invariant 
mass distribution has been taken into account. Changing the background de- 
scription, the degree of polynomial, the top mass peak shape from a Gaussian 
to a Breit-Wigner or to the sum of two Gaussians with common means lead to 
a systematic error on the top mass of 90 MeV. 

The effects of initial and final state radiation have been computed in the 
same way as before. The systematic error due to ISR is found to be negligible 
and the error due to FSR is found to be 200 MeV. 

Since the light quark jets are calibrated in-situ, a negligible systematic error 
in the top mass measurement results from the uncertainty in the light quark jet 
energy scale. For the b-jet energy scale, a 1% miscalibration induces a top mass 
shift of 700 MeV, as was discussed in section 12.3.21 for the kinematic fit with 
fixed cone size. 

A study was performed to investigate whether one could reduce the b-jet 
energy scale systematic error by calibrating the b-jets using the same in-situ 
calibration obtained for the light quark jets. Proceeding this way will increase 
the statistical error, due to the data based calibration, but also introduce a 
systematic shift of the top mass due to any differences of energy losses between 
light quark jets and b-jets. These differences are expected due to physics effects, 
and also due to detector and reconstruction effects. The systematic error due to 
the b-quark fragmentation parameter et was estimated as before and was found 
to be 50 MeV. The b-jet energy scale depends also on the branching fraction 
of the semi-leptonic b-hadron decays, due to the presence of the undetected 
neutrino. To estimate the influence of the imprecise knowledge of the semi- 
leptonic decay fraction of b-hadrons (which is known with an accuracy of 7%) , all 
semi-leptonic branching ratios were scaled by 1.07 in the Monte Carlo, resulting 
in a top mass shift of 60 MeV. Due to the fixed jet cone reconstruction procedure. 
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Fig. 18. Reconstructed top mass. The full fine represents the result of a fit by a 
Gaussian plus a 4th order polynomial. 



the b-jet energy scale depends on the parton shower evolution and in particular 
on the as evolution parameters. Introducing the b-quark mass corrections into 
the default as{pT) evolutions in showers led to a systematic error of 90 MeV. 
Combining these uncertainties in quadrature would give a total systematic error 
on the top mass due to the differences of the physics effects between light quark 
jets and b-jets of 130 MeV. While this result is very encouraging, it must be 
stressed that these comparisons were made using a fast, parameterized Monte 
Carlo description of the ATLAS detector. Further study will be performed with 
full GEANT-based simulation of the detector to increase the confidence of this 
potential reduction in the systematic error. 

The shape of the W signal used to rescale the b-jet energy is another source 
of systematic uncertainty. The systematic error was estimated by taking into 
account the asymmetric shape of the W signal distribution, changing the back- 
ground shape under the W peak and the fitting region. These sources account 
for a systematic error on the top mass of 100 MeV. 

The various contributions to the top mass systematic error for the contin- 
uous jet analysis arc summarized in table |^ Adding the various contributions 
in quadrature leads to a total systematic error on the top quark mass of the 
order of 1 GeV, dominated by the uncertainty in the external calibration for the 
b-jets. Should it be possible to realize the improvement suggested by the study 
of b-jet calibration using the in-situ light quark calibration, the total error could 
be reduced to of order 400 MeV. 



2.5 Summary 

Three methods to determine the top quark mass in the lepton plus jets channel 
have been presented in chapters 2.2, 2.3 and 2.4. In the first method (2.2) the top 
mass is extracted from the invariant three jets mass of the hadronic top decay. 
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Source 


Smt (GeV) 


Range of jet cone sizes 


0.25 


dependence 


0.2 


Signal and background shape 


0.1 


ISR and FSR 


0.2 


External b-jet calibration 1% 


0.7 


Internal b-jet calibration 




Physics effects 


0.13 


W signal shape 


0.1 



Table 6. Summary of the systematic errors for the top mass in the inclusive 
lepton plus jets sample, when the top mass is reconstructed using a continuous 
jet definition. For more details, see the text. 



in the second method (2.3) the top mass is determined from a kinematical fit 
of the entire ti decay, and in the third method the top mass is determined 
from a kinematic fit and using a continuous jet algorithm. The main sources of 
uncertainties arise from final state radiation and b-quark jet energy scale. It was 
shown that the contribution from light quark jets energy scale to the systematics 
errors can be reduced to a negligible level using an in situ calibration. Provided 
that the b-quark jet absolute energy scale can determined within 1%, the top 
quark mass can be measured with a precision at the level of 1 GeV in one year 
of LHC running at low luminosity (per 10/6^^). 

2.6 Top mass measurement using large px events 

In this section, we present an alternative method which uses a special sub-sample 
of the single lepton plus jet events where the top has high transverse momentum, 
for example px > 200 GeV. 

In this topology, the two quarks are produced back-to-back, and the daugh- 
ters from the two top decays would appear in distinct hemispheres of the de- 
tector: the "hadronic" one from the decay t Wb — + jjb, and the "leptonic" 
one from the decay t Wb — > Ivb. Due to the high pT of the event, detector 
systematics as well as backgrounds from other processes are expected to be very 
small. The distinct feature of these events which is exploited here, is the fact that 
due to the high px the three jets from the hadronic top decay tend to overlap 
in space, as shown in figures [T9I2UI and 12 II 

Therefore one could reconstruct the top mass without using the jets as in the 
methods described in the previous section, but from summing up the individual 
calorimeter towers over a large cone {AR in [0.8-1.8]) around the top direction. 
The top direction itself can be determined in two ways: i) as opposite to the top 
direction reconstructed in the leptonic decay, where the missing energy in the 
event is used to reconstruct the neutrino, and ii) as the direction of the invariant 
mass of the three jets in the hadronic top decay. Figure 1^ shows the percentage 
of the generated events with all the three jets from the hadronic top decay lying 
within a distance AR from the top quark direction. 
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Fig. 19. AR distance between the top Fig. 20. AR distance between the top 
quark and the furthest quark from the quark and the closest quark from the 
hadronic W decay. hadronic W decay. 




1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 



' AR(&h -b-jetsj 

Fig. 22. Percentage of the events hav- 
Fig. 21. AR distance between the top ing all the three jets from the hadronic 
and the b-quark at parton level. top decay within a cone aperture AR 

from the top, at parton level. 
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Using directly the calorimeter towers, avoids problems with the jet recon- 
struction and energy calibration, one of the major source of systematic errors in 
the top mass measurement, and introduces a different set of systematic errors. 
This method is therefore very interesting and useful for the final combined top 
mass determination by ATLAS. 

2.6.1 Event selection and reconstruction Samples of high pT tt events 
were generated with a cut of 200 GeV in the center of mass of the hard scattering. 
The cross-section of this topology corresponds to about 2% of the total ti cross- 
section. Events are selected to pass the trigger selection and by requiring one 
isolated lepton with px > 30 GeV and I77I < 2.5, the transverse missing energy 
greater than 30 GeV (i?™'** > 30 GeV), and at least four jets reconstructed 
using a cone of AR = 0.4 with pt > 40 GeV and \ri\ < 2.5, of which two must 
be tagged as b-jets. 

The overall efficiency is 9%, resulting to ~15000 selected events per (10 
fb~^) ^. Due to the high px of the event and the requirement for two tagged 
b-jets, the background (mainly W-|-jet, WW or QCD events) is reduced into 
negligible levels and therefore not discussed further. 

For the events passing the preselection cuts described above, the top quark 
direction was determined as described. First the hadronic W invariant mass was 
reconstructed from the two highest pr non b-tagged jets. Combinations where 
rrijj = My/ ± 20 GeV were selected. The two jets were combined with the closest 
b-jet to reconstruct the top. Finally, the reconstructed top pt was required to 
be above 235 GeV. After all the cuts ^3600 events remain per 10~^ fb, with an 
overall efficiency of 2%. 

Once the top direction is determined, the invariant mass of all the calorimeter 
towers {Arj x A<j) = 0.1 x 0.1) around this direction is evaluated according to the 
formula: 

n(AR) n(AR) 
i=l i=l 

where Ei is the total calorimeter energy in the i-th tower evaluated in elec- 
tromagnetic scale, and Pi is its three momentum vector. The index n{AR) runs 
over all the towers within the selected cone radius. This invariant mass is directly 
proportional to the top quark mass: rriciust = '^ctesr 

Figure 1231 shows the reconstructed fn^iust invariant mass for a cone size of 
Z\i?=1.3. A clear Gaussian distribution is observed with the peak value around 
the nominal top mass. Fitting the peak region with a Gaussian, we obtain a peak 
width of 9.6 GeV, comparable to that obtained with the jet method. 

As shown in figure !^ more than 80% of the events where the three jets are 
at AR < 1.3 from the top quark direction are selected. 

The invariant mass w*"^^^, is evaluated using various cone sizes in the range 
from AR— 0.8 to 1.8, as there is no reason a priori to select a given value. On 
the contrary, the resulting invariant mass has to be independent from the cone 

^To save computing time, the studies presented here are done using the muon channel only, 
and assuming similar efficiencies for the electrons. Several ATLAS studies have demonstrated 
that this is a good approximation when working with high px objects as in this analysis. 
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size used. For each AR cone size, a Gaussian is fitted around the peak position 
in the invariant mass distribution. The fitted peak values for the different 

cone sizes are shown in figure YI^ The variation observed can be attributed to 
the Underlying Event contribution which is added for each calorimeter tower, re- 
sulting in an increased invariant mass value as the cone size increases. A method 
to evaluate the UE contribution in each tower follows. 

2.6.2 The Underlying Event (UE) estimation The Underlying Event 
contribution per calorimeter tower (UEtower) was estimated from the same high 
Pt top sample. It represents the average transverse energy Et deposited per 
calorimeter tower in each event, once all the towers related with the high pT 
products are excluded. The values as well as the number of towers used in each 
case have been computed for different rapidity regions j22j. An average over all 
rapidity and isolation cut range, gives a value of U Etower—4=4:7 .5 MeV, which is 
subtracted from the energy of each tower in the calculation of rriciust- 

In figure 1241 the invariant mass is shown after the the UEtower is 

subtracted. The resulting values are now independent of the cone size, with 
an average value of 159 GeV and with all values within ±0.15%. Varying the 
UEtower by ±10%, a cone size dependence is again observed, raising to a bit less 
than ±2%, which demonstrates that the value used is the correct one and gives 
the precision required for the target top mass measurement error. In ATLAS, 
once the real data become available, the UEtower will be calculated in situ as 
done here, but also using other event samples, resulting to an overall error of 
about 10%. 

2.6.3 Mass scale calibration After the UEtower contribution is subtracted, 
the reconstructed rn^^^^^ invariant mass has become independent of the cone size, 
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but the resulting values are now ~9% lower than the generated top quark mass. 
This a priori is expected as no particular mass scale or absolute energy scale has 
been used so far. As a first attempt the w*"^^,^ values could be calibrated using 
the Monte Carlo data. However doing so the method will be dependent on the 
exact modelling of the process and won't be anymore a "direct" measurement of 
the top mass. The best way is to obtain a mass scale calibration using the data 
themselves. 

The method studied here is to apply the same reconstruction procedure, with 
the same UEtower, to other known particles with well measured masses and 
extract from there the necessary mass scale calibration factors. In our case the 
easiest way is to use the inclusive top sample and apply the same reconstruction 
method to the W mass reconstruction. This sample offers high statistics, and 
has practically the same event topology as the high pT sample. Rescaling the 
corresponding rn^^^^ invariant mass values to the nominal mass of the W, we 
can obtain the mass scale calibration coefficient Ctop averaging all the cone sizes 
AR. Finally, to determine the top mass, the average value of all the rn°j^^^ values 
after calibration is used. Fixing the mass scale with the W, and then transferring 
the results to the top, it implies that the same calibration is used for the both 
the light quarks and the b-jets. 

Events were generated in the single lepton plus jets topology, without pt cut 
applied to the hard scattering process. The large statistics available, allow to 
apply tight cuts in order to select events where the two jets from the hadronic 
W decay are close in space (as in the high px sample) and at the same time far 
away from the b-jet. Events were selected by requiring an isolated lepton with 
PT > 20 GeV and I??! <2.5, EJJ^^"'' > 20 GeV, at least four jets (reconstructed 
in a cone of AR—O.A), with pT > 40 GeV and \ri\ <2.5, of which two are tagged 
as b-jets. In addition, the distance between the two highest non b-tagged jets, 
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Fig. 27. Reconstructed spec- 
trum obtained using a cone size of 
AR=1.3 around the top direction. 



Fig. 28. The ra°j^^^ invariant mass de- 
pendence on the cone size before and af- 
ter the UEtower is subtracted. 



should be AR <1.3 and the two b-jets of the event should be at a distance 
AR >2.0 away from the reconstructed W direction. 

The two highest pt non b-tagged jets were used to reconstruct the W and 
find its direction. Then the m^^^ invariant mass was calculated around this 
direction subtracting from each tower the same UEtower as calculated before. 
In figure 123 the fitted value of the invariant mass rri^^^^ are shown. As for the 
top, the reconstructed values after the UE is subtracted become independent 
from the cone size (within ±0.7%) and about 7.7 GeV lower than the nominal 
W mass. 

Figure shows the resulting ?ti*°^^j values after the mass scale calibra- 
tion is applied. The variation between the points is ±0.22%. As an example, 
for Z\i?= 1.2 (1.3) the nri°i^^gt invariant mass after the UEtower subtraction was 
158.9(159.0) GeV, and after applying the calibration becomes 175.7(175.8) GeV. 
Taking the average of the calibrated nri°j^^^ invariant mass for all cone sizes, a 
value for mt — 175.9 GeV is obtained, which is within 0.5% from the generated 
top quark mass. 

2.6.4 Full simulation results The results presented so far were obtained 
using the fast detector simulation j^. The same analysis was repeated with a 
sample of full (GEANT-based) simulated events of the ATLAS detector. Figure 
123 shows the reconstructed to*;^^^ invariant mass spectrum for a cone size of 
AR=l.'i. The variation of the fitted values for different cone sizes is shown in 
figure |2i 

Although the peak values in this case are lower than those of the fast simula- 
tion, the overall variation for the same cone size range stays about the same. The 
difference between the fast and full simulation can be attributed to the shower 
shape development which is not included in the fast simulation. 
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The UEtower was evaluated following the same procedure as before. The 
average UEtower is now 42.5 MeV much lower than the fast simulation. 
Since now there are more calorimeter towers contributing but with lower energy 
in each, compared to the fast simulation case where the energy for each particle 
is deposited to a single tower. This value was used for all the full simulation 
studies described below. 

The rn^^st invariant mass after the UEtower subtraction is shown in figure 
1281 as a function of cone size used. As expected, it remains basically independent 
from the cone size, but lower by 24.6% from the generated top mass. The vari- 
ation observed, ±0.9%, is bigger than with the fast simulation sample, and can 
be attributed to the poor quality of the fits due to the lack of statistics. Using 
only the points up to Z\_R=1.4 the average value is 133 GeV, with a variation of 
±0.2%. 

The same mass scale calibration procedure was used, with a sample of 30000 
fully simulated inclusive tt events, applying the same cuts and reconstruction 
procedure. The fitted m^^^^ peak values for different values of AR are shown in 
figure !^ before and after the UEtower is subtracted. After the UE contribution 
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is subtracted, the resulting rn^^^^ invariant mass remains independent from the 
cone size but 23.5% below the nominal W mass. The mass scale is determined as 
before, and "rn^^^t invariant mass values after the calibration factors are applied 
are shown in figure 1^ The resulting values are constant within ±1.4% from the 
generated top quark mass, for the whole range of cone values used. The average 
rut is 172.6 GeV, which is 1.4% below the generated top mass, and with all points 
within ±0.9%. Restricting to the range up to AR <1.4, where basically we run 
out of statistics, the mt changes to 173.8 GeV, with the points now having a 
spread of ±0.3%. 
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In summary, the full and fast simulation data show the similar results and 
confirm the competitiveness of the proposed reconstruction method. However, 
further studies with larger statistics samples should be made in particular for 
high AR values. 

2.6.5 Systematic uncertainties Several studies have been performed which 
cover most of the possible systematic errors. Studies requiring large statistics and 
several settings of generator parameters were performed with the fast detector 
simulation, while for the cases where the exact detector response is important 
samples of full detector simulation were used. In some cases where large com- 
puting effort was required only the sensitivity of the results was investigated 
without going to details. More information on these studies can be found in |22|. 

To study the linearity of the reconstructed top mass several samples of high 
Pt top events with different input top quark mass in the generator from 160 
GeV to 190 GeV were produced and analyzed in exactly the same way. For all 
the samples, the same UEtower and mass scale calibration factors obtained as 
explained before were used. In the mass range 170-180 GeV, the reconstructed 
top mass is in very good agreement with the generated value. For larger values 
a bigger error is observed which is somehow expected as the event environment 
changes and the UE and mass scale calibrations are not optimal anymore. 

The sensitivity to initial (ISR) and final (FSR) state radiation was studied in 
the same way as in other analysis. Samples of high top events were generated 
with the ISR or FSR contributions switched off at the generator level, and the 
analysis was repeated keeping the same U Etower estimate and the mass scale 
calibration factors. Doing so, and for a cone radius AR — 1.3 a shift in the 
reconstructed top mass of 0.7(0.3) GeV is observed when ISR(FSR) was not 
present. It has to be pointed out this error is a very pessimistic approach, as the 
exact level of the ISR(FSR) contributions in the events will be measured and 
known at LHC to about 10% level and the generators will be correctly tuned to 
this. Therefore as in the other analysis, the final error quoted for the top mass 
is 20% of the total mass shift observed, equal to 0.1 GeV. 

The sensitivity on the b-quark fragmentation was studied by generating sam- 
ples where the parameter in the Peterson formula was varied within its error 
currently at 0.0025 The top mass was reconstructed in each sample using 
the same U Etower and mass scale. The observed top mass shift among the sam- 
ples is quoted as the error due to this effect. As an example, for a cone size 
of Z\i?=1.3 the mass shift was 0.3 GeV. Similar values obtained for other cone 
sizes. 

The UE energy estimate per tower plays an important role in this top re- 
construction method. To evaluate the sensitivity of the reconstructed mass due 
to this, the value calculated (447.5 MeV for the fast and 42.5 MeV for the full 
simulation data) was varied by ±10% and the top reconstruction and the mass 
scale calibration was repeated each time. As shown in figure|^the reconstructed 
values stay well within ±1% of the generated value. 

To study the contribution of a possible calorimeter mis-calibration in the top 
mass measurement, the energy in each tower was varied according to a Gaussian 
with different values of sigma from 1 to 5%, well beyond the expected reach by 
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Fig. 31. Reconstructed m°i^^^ mass de- 
pendence on the UEtower energy esti- 
mate. At each case the mass scaling fac- 
tors have been recalculated as explained 
in the text. 



Fig. 32. Reconstructed m°if^^f mass val- 
ues before and after the rescaling from 
the W reconstruction in the inclusive 
sample for three values of e/h. 



ATLAS. The analysis was repeated in each case keeping the UEtower unchanged. 
For a cone size of AR = 1.3 the observed shift in the reconstructed top mass 
is 0.6(1.2) GeV when a mis-calibration of 1(5)% is added (0.7(1.3) GeV in full 
simulation). This value however is rather pessimistic, since whatever the cell 
mis-calibration would be the UEtower would be computed accordingly, and the 
whole error will be incorporated into the mass scale calibration. 

Tests with the ATLAS calorimeter prototypes have shown that the combined 
e/h is different than 1.0 23- Moreover it was demonstrated that in both fast and 
full detector simulation programs the e/h effect on the calorimeters is not cor- 
rectly treated . To study the sensitivity of the reconstructed top mass due to 
this, given that in this method the individual calorimeter towers that combine 
information from several calorimeters is used, several fast simulation samples 
were generated where the energy deposited by hadrons has been corrected ac- 
cording to the Groom's formula |2S] thus simulating values of e/h between 1.0 
and 1.63. The whole analysis was repeated in each case. Changing e/h to 1.63, 
a rather extreme value, the UEtower estimate changes from 447.5 MeV to 417.5 
MeV, and the reconstructed top mass by 0.7%, as expected from the mass scale 
calibration procedure. 

Since for the top mass reconstruction the whole calorimeter volume is used, 
either to calculate the UE contribution either to evaluate the top mass, the 
presence of electronics noise may have an impact on the results. To study this, a 
test using the full simulated data was performed, where the electronics noise was 
added as Gaussian noise to the cell energies according to the expected values for 
each calorimeter. Repeating the same procedure as before, the UEtower value 
changes from 42.5 MeV to 16 MeV, a surprising low value. Following the same 
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procedure and applying the mass scale calibration, the average reconstructed 
top mass, is 175.17 GeV but shows a variation within ± 1.2% for the full range 
of cone sizes used between 0.9 and 1.6, partially enhanced due to the lack of 
statistics. Further studies are needed to fully understand the effect, once the 
detector is built. 

2.6.6 Summary The method presented here uses a special sub-sample of the 
single lepton plus jet events where the top has high transverse momentum. The 
case with pT > 200 GeV was studied here, which due to the high tt production 
at LHC offers good statistics per year at low luminosity (lOfb^^). 

After some pre-selection cuts aiming to efficiently select the hadronic top de- 
cay products with minimal contribution from the background, the top direction 
is found using jets similarly to the inclusive sample. Then, using the unique fea- 
ture of these events, where the hadronic top decay products are well collimated 
in space, all the calorimeter towers within a cone of radius ARdust around the 
reconstructed top direction are summed, forming a top invariant mass. It was 
then shown that the reconstructed invariant mass becomes independent of the 
cone radius used (varied between 0.8 and 1.8) once the underlying event con- 
tribution is subtracted. Finally, the mass scale was determined by applying the 
same reconstruction method for the W jj decay in the top events of the inclu- 
sive sample. After recalibration, the results obtained with both the fast and full 
detector simulation are comparable with the other methods, and demonstrate 
that the top mass can be reconstructed with an accuracy of 1%. 

Using this method the reconstructed top mass is sensitive to a different set 
of systematic errors, a first study of which as was performed. From the results 
obtained, the major contribution is in the mass scale calibration procedure, with 
all the errors stay below the 1% level. Table [7| summarizes results. For some of 
the studies described above, the final systematic error once the detector and the 
data are available will be incorporated in the mass scale calibration procedure, 
therefore are not included as individual lines in the table. 

\Amt\ (GeV) Smt (GeV)~ 



Initial state radiation 0.7 0.1 

Final state radiation 0.3 0.1 

b-quark fragmentation 0.3 0.3 

UE estimate (±10%) 1.3 1.3 

mass scale calibration 0.9 0.9 



Table 7. Top mass shift (\Amt\) and quoted systematic error on mt (6mt) due 
to various sources of systematic uncertainties. 



3 Top mass measurement in the dilepton channel 

The dilepton events can provide an indirect measurement of the top quark mass. 
The difficulty comes from the fact that, in principle, one cannot fully reconstruct 
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the top decays due to the presence of undetected neutrinos in the final state. 
For the determination of the top mass, previous methods have exploited the 
correlation between the top mass and kinematic quantities, such as the mass 
mih of the lepton-b-jet system 0. 

Here, in a first step, assuming a value for the top mass, it is proposed to 
reconstruct the top decays by solving the set of equations describing the kine- 
matic constraints of the decays. Then, to determine the top mass, the solutions 
obtained for different input top masses will be compared to the data |26| . 

3.1 Event selection 

The dilepton events are characterized by two high px isolated leptons, large 
traverse missing energy i?™'"^ and two jets coming from the fragmentation of 
the b-quarks. Taking into account the branching ratio, about 400000 dilepton 
events can be expected for integrated luminosity 10/6^^. 

The background is coming mainly from Drell-Yan processes and Z tt 
associated with jets, and from WW+jets and bb production. Events are selected 
by requiring two opposite sign isolated leptons with px > 35 GeV and pr > 25 
GeV respectively and |r?| < 2.5, > 40 GeV, and two jets with > 25 GeV. 

After event selection, 80000 signal events are left, with a signal over background 
ratio around 10. 



3.2 Method for the final state reconstruction 

For the determination of the momenta of both the neutrino {f) and anti- neutrino 
(p), it is assumed that the masses of the top and anti-top are known. The re- 
construction algorithm is based on solving a set of equations coming from the 
kinematic properties of the conservation of momentum and energy |^ . The set 
of equations consists of six equations for six unknown components of momenta of 
neutrino and antineutrino. First two equations describe conservation of transver- 
sal momentum of the ti system, assuming that this momentum is 0. The other 
equations constrain invariant masses of both lepton-|-neutrino systems to the 
masses of W'^ and W~ bosons, and masses of both lepton-|-neutrino-|-jet system 
to the masses of top and antitop quarks. All W^, , top and antitop masses 
are assumed to be known. 

After some derivations |2|^, the following two linear equations with the un- 
knowns , Py, and , are obtained: 



' St 
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+2 (^pi -Ks + p'',+kl + kl + kl + Mt + El - kl 

-pi -P'y -pI -p\ +2 I "^^^ \^ " +k,E, + k,p'^ + k2p'y\ 

-2 + M + P\p\ + Mv^+) = 

Where: Ep is the energy of particle p, Mp is a function of the mass of particle 
p, ki is a function of momenta of leptons and b-quarks, p^ represents the i-th 
component of momentum of particle p, I represents either e+ or I represents 
either e~ or u is either Vf, or i/^, j7 is either 17^ or P^. 

Additional derivations lead to one quartic equation with only one unknown, 
which is analytically solved. All the derivations were performed using a software 
tool for symbolic algebraic manipulations. 

The remaining components of both neutrino and anti-neutrino momenta can 
be easily computed. Finally, the complete kinematic reconstruction can be per- 
formed. 

The reconstruction algorithm can provide no solution or more than one solu- 
tion. In the first case, the right-handed sides of the two equations from the initial 
set of six equations describing momenta conservation are varied in the range [- 
250 GeV:+250 GeV] starting with until an acceptable solution is found. The 
solubility of the system is improved from 88% to 97.6% (this means that the tt 
decay is reconstructed for 97.6% of the events). In the second case, the choice 
of the solution is based on the computing of weights for known distributions of 
various kinematic quantities of the ti decay "26' . The right solution is chosen in 
73% of the events. 

Therefore, the reconstruction algorithm exhibits an efficiency of 97.6% with 
a purity of 73%. 

3.3 Top mass determination 

It has just been demonstrated that the entire tt decay can be reconstructed by 
assuming a value for the top quark mass. For the top mass determination, the re- 
construction algorithm will be fed with various top masses and the corresponding 
solutions will be compared to the data. 

3.3.1 Method The method was tested using samples of events containing 
approximately the same amount of events that will be collected during one year 
of running at low luminosity (per 10/6~^), after selection cuts are applied. 

The principle of the determination of the top mass is the following: for each 
event, one tries to solve the equations for various input top masses, and to 
compute the weight of the best solution for a given top mass value. If the input 
top mass value is quite different from the correct value, no solution may be 
found, or the solution will have a small weight. 
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Fig. 33. Mean weight as a function 
of the input top mass. The maximum 
mean value gives the top mass. 



Fig. 34. Mean weight as a function of 
the top mass, with ISR switched off. 



For each input top mass value, a mean weight over the entire set of events 
is computed. The top mass is given by the value having the maximum mean 
weight. 

The weight mean value as a function of the input top mass is represented on 
figure ISSl As expected, the curve is peaked around the generated top mass value 
of 175 GeV. The maximum mean weight is obtained by fitting this curve with a 
quadratic function, leading to a reconstructed top mass in agreement with the 
generated value. The error on the reconstructed value is 0.3 GeV. This error 
takes into account both statistical effects and effects due to the reconstruction 
method itself. 



3.3.2 Systematic uncertainties The effect of the systematic uncertainties 
sources on the top mass determination has been checked following the methods 
described in section 2. For initial and final radiation, top mass shifts were de- 
termined by switching ISR and FSR off separately in the Monte Carlo, and the 
corresponding error was obtained by taking 20% of the mass shifts. The b-quark 
jets calibration was assumed to be known within 1 %. For the b-quark fragmen- 
tation parameter, a mass shift was computed between the top mass obtained 
with the default value Sb — 0.006 and with Sb = 0.0035. The error due to the 
parton distribution function was estimated by measuring the top mass shift be- 
tween a sample of events simulated with the default set and a sample of events 
simulated with another set. The values are summarized in table |S| 

An example, for ISR switched off, of the mean weight as a function of the 
input top mass is shown on figure 1341 The values of the estimated systematic 
errors listed in table |H| are not large. This is due to a very positive aspect of the 
reconstruction method. For a given systematic uncertainty source, the curve of 
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source of uncertainty 



\Amt\ (GcV) 5mt (GeV) 



Statistics and reconstruction method 

b-jet energy scale 

b-quark fragmentation 

Initial state radiation 

Final state radiation 

Parton distribution function 



0.6 0.6 

0.7 0.7 

0.4 0.1 

2.7 0.6 

1.2 1.2 



0.3 



Table 8. Top mass shift Anit and resulting systematic error on mt{Smt) due 
to the various source of systematic errors, in the dilepton channel. 



the mean weights versus the input top mass is modified in two ways compared to 
the curve obtained with the default sample. The mean weights are smaller, giving 
a maximum mean value smaller than the initial value, and the peak is shifted 
giving a corresponding shifted top mass. This second effect only is relevant in 
the systematic uncertainty studies. 

3.4 Summary 

It was shown that, assuming a mass for the top quark, the final state topology of 
dilepton events can be fully reconstructed by solving a set of equations describing 
the kinematic constraints of the tt decay. The decay reconstruction algorithm 
has high efficiency and purity. A step further, for the determination of the top 
mass, consists in feeding the reconstruction algorithm with different input top 
masses and to compare the solutions with the data. 

There is also a possibility to consider more than two jets in final state. In this 
case one has to solve the set of equations for all 2-jets combinations. Surprisingly, 
this has also no impact on the estimation of the top mass value, however, it is 
one of the subjects to be studied yet. 

A preliminary study of the systematics uncertainties shows that the top mass 
can be extracted with a reasonable accuracy, at the same level as other tech- 
niques. This method can therefore provide a useful input for the combined AT- 
LAS top mass measurement. 

4 Top mass measurement in the six jets channel 

The all jets channel final state topology consists, in the absence of initial or 
final state radiation, of six jets (including two 6-jets), no high px leptons, and 
small transverse missing energy Et- With no energetic neutrinos in the final 
state, the all hadronic mode is the most kinematic-ally constrained of all the tt 
topologies, but it is also the most challenging to measure due to the large QCD 
multijet background. Nevertheless, at the Fermilab Tevatron Collider both the 
CDF and D0 collaborations have shown that it is possible to isolate a tt signal 
in this channel [271 05| . The CDF collaboration obtained a signal significance 
over background of better than three standard deviations '27' by applying simple 
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selection cuts and relying on high 6-tagging efficiency. To compensate for the less 
efficient 6-tagging, the D0 collaboration developed a more sophisticated event 
selection technique based on a neural network j28j . 

The potential of the ATLAS detector to study the all hadronic decays of tt 
pairs has been explored. In the search for an optimal strategy for signal extrac- 
tion from background, the kinematic properties of both signal and background 
events are investigated, and a kinematic fit of selected events is performed. Fi- 
nally, a clean sample is obtained by selecting events in which both reconstructed 
top and antitop quarks have a high traverse momentum [px > 200 GeV). This 
subsample is then used for the reconstruction of the top mass. Nevertheless, the 
top mass is reconstructed in the inclusive sample as well. 

4.1 Signal selection 

Taking into account the branching ratio, the next-to-leading order cross-section 
prediction for the all jets channel is 370 pb. Therefore, for an integrated lu- 
minosity of 10 fb~^, one can expect 3.7 million ti pairs with this final state 
topology. 

The main source of background is QCD multijet events, which arise from 
2^2 parton processes (g^Qj -> q^qJ, q^g q^g, gg gg, qq -> g_g, qiqi -> q^qj, 
gg — > qq) convoluted with parton showers. The heavy-flavor (cc, bh, tt) content in 
a QCD multijet sample stems from direct production (e.g. (j^gi QjQj^ 99 ^ 
gluon-splitting (where a final state gluon branches into a heavy quark pair), 
and flavor excitation (initial state gluon splitting). In the analysis that follows, 
ti production was excluded from the QCD background processes. The QCD 
background was generated with a pT cut on the hard scattering process above 
100 GeV, resulting in a production cross-section of 1.73 fib. Processes involving 
the production of W and Z bosons (with their subsequent decay into jets) were 
not included since their contributions are small compared to the QCD multijet 
background. 

As the first step in the selection of the all hadronic tt topology, events were 
required to have six or more reconstructed jets, of which at least two must 
be tagged as 6-jets. Jets were reconstructed using a fixed cone algorithm with 
AR=OA. Jets were required to have pt greater than 40 GeV, and to satisfy jryj < 
3 (It^I < 2.5 for 6-jet candidates). The efficiencies for these selection criteria for 
both tt signal and QCD multijet background are 2.7 % and 0.011 % respectively, 
resulting in a signal over QCD background of 1/19, indicating that these simple 
selection cuts can already reduce the multijet background to manageable level. 

4.2 Signal and background kinematic properties 

Further progress in enhancing the S/B ratio could be sought using variables that 
provide discrimination between the signal and the QCD background. Therefore, 
some kinematic variables sensitive to the energy flow in the event, additional ra- 
diation and event shape (including several variables used in the neural network 
analysis of the D0 collaboration |5S]) were examined. Those variables include: 
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Fig. 35. Distributions for tt signal (hatched) and QCD multijet background (un- 
hatched). Left: the Ht, -ff^"''^*'*, E'^-''^** and y/I distributions. Right: apla- 
narity, sphericity, centrality and AR; 
details. 
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Ht'- the sum of all jet transverse energies in the event ( ^ Et ). 

i=i 

: Ht without the transverse energy of the two leading jets. 
ji*^jets 1,2. ^j^g transverse energy of the two leading 6-jets. 

the invariant mass of the jets in the final state. 
A: the aplanarity, fQi, calculated from the normalized momentum tensor. 
S: the sphericity, §((9i + Q2), calculated from the normalized momentum tensor. 

C: the centrality, Ht/He, where He — J2 -^j ^^e sum of all the jet total 

energies. The centrality characterizes the transverse energy flow. 
ARjj™: the minimal separation between two jets in r/-0 space. 

The first four of these variables are related to the energy deposition in the 
event, while the others are more related to the event shape or topology. The 
normalized distributions for these variables, for p:^^* > 40 GeV, are plotted for 
tt signal and QCD background in figure (left plot for the first four variables 
and right plot for the others). It can be seen that the variables sensitive to 
the event shape provide a somewhat better discrimination between the signal 
and background. However, it is clear that none of these variables provides at 
the LHC the clear discrimination which was observed at the Tevatron energy 
|28| . Therefore, it would appear difficult to select a relatively clean signal based 
on cuts on these variables, or even the use of a more sensitive cut based on a 
multivariate discriminant, where the variables are treated collectively |28l 129) . 
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4.3 Final state reconstruction with a kinematic fit 

The key feature distinguishing top quark events from QCD multijct background 
is the fitted mass obtained from the least-squares kinematic fit of the events to 
the tt decay hypothesis p!!!. In order to simplify the analysis, massless jets have 
been assumed and the error on the measured jet direction was neglected with 
respect to the error on the measured jet energy. 

The reconstruction algorithm and fitting procedure proceed in two steps. 
First, the two W jj decays are reconstructed by selecting di-jet combinations 
from jets not tagged as 6-jets. This is done by minimizing a Xw function 

Next, the two W jj candidates are combined with the 6-tagged jets to form 
the top and antitop quark candidates {jjb combination). The energies of the b- 
and 6-jets are constrained by minimizing a Xt function |31| . There are two ways 
to associate the b-tagged jets to the reconstructed W bosons. The association 
giving the smallest value of Xt chosen. After the event reconstruction and 
fitting procedure, additional qualitative cuts are applied |31j . 





after selection 


after kinematic fit 


within the window 




cuts 


and x^ cuts 


130-200 GeV 


tt{%) 


2.7 


0.3 


0.18 


QCD (%) 


0.011 


0.00017 


0.000007 


S/B 


1/19 


1/2.6 


^6/1 



Table 9. Efficiency for tt signal and QCD multijct background after applying 
various level of cuts, and for pl^* threshold > 40 GcV. The last row shows the 
resulting signal to background ratio. 

Table |51 presents the efficiency and S/B ratio for tt signal and QCD multijet 
background after selection cuts are applied, after the kinematic fit procedure, 
and after the additional requirement that the reconstructed top and antitop 
quarks masses lie within the window 130-200 GeV. The kinematic fit and limits 
on the top (antitop) mass significantly improve the value of S/B ratio so that, for 
top masses within the 130-200 GeV window, the S/B ratio is « 6. This result is 
obtained with a large statistical error on the remaining background (40%), and 
a more accurate determination would require generation of significantly larger 
Monte Carlo background samples. 

4.4 High transverse momentum tt events 

The signal over QCD background can be further improved by restricting the 
analysis to a sample of high transverse momentum tt events where both recon- 
structed top and anti-top quarks have pr > 200 GeV. To study this sample, tt 
signal and QCD background events were generated with a px cut on the hard 
scattering process above 200 GeV. The corresponding cross-sections are 53.5 pb 
for signal and 86.1 for QCD multijet background. 
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Fig. 36. Invariant mass distribution of the accepted jjb combinations for the 
high pt sample, normalized to an integrated luminosity of 10 fb~l. The shaded 
area shows the QCD multijet background. 



4.4.1 Top mass reconstruction For the highpx sample, the same selection 
cuts were applied as for the inclusive sample. After the kinematic fit and the 
requirement that both reconstructed top and anti-top quarks have px > 200 
GeV, the selection efficiencies and S/B ratios are given in table 1101 

The invariant mass distribution of the accepted jjb combinations and for the 
QCD background (the shaded area) is shown in figure Within the window 
130-200 GeV the signal over background ratio is « 18. 

The distribution fitted by a Gaussian leads to a reconstructed top mass 
consistent with the generated value with a peak width of 10.1 GeV. For an 
integrated luminosity of 10 fb a sample of 3300 events would be collected with 
fully reconstructed top and antitop quarks with pt > 200 GeV. This number of 
events would lead to a statistical error of (5mt(stat)=±0.18 GeV. It can be noted 
that this clean sample could be used for the study of differential distributions 
for both top and anti-top quarks |31| . 
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130-200 GeV 


eu{%) 


0.68 


0.63 


eQCD{%) 


0.00041 


0.000021 


S/B 


1/1 


« 18/1 



Table 10. Efficiency for high pT tt signal and QCD multijet background where 
reconstructed top and antitop quark both have px > 200 Ge V, for different cuts 
applied and for p^^* threshold > 40 Ge V. The last row shows the resulting signal 
to background ratio. 
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4.4.2 Systematic uncertainties The systematic uncertainties have been 
treated in a similar way as in the inclusive lepton plus jets channels. It was 
assumed that jet energy scale for both light quark and b-quark jets will be 
known at the level of 1%. For the b-quark fragmentation parameter Sb, a top 
mass shift was determined between the top mass obtained with the default pa- 
rameter {sb = 0.006) and with Eb = 0.0035. For initial and final state radiation, 
mass shifts were obtained between ISR and FSR switched on and off separately. 
The resulting systematic error was taken by considering 20% of the mass shifts. 
The results are summarized in table 1111 for the high px sample. 



Systematics 


5mt (GeV) 


Light jet energy scale 


0.8 


b-jet energy scale 


0.7 


b-quark fragmentation 


0.3 


Initial state radiation 


0.4 


Final state radiation 


2.8 



Table 11. Systematic error on mt{Smt) due to the various source of systematic 
errors for the high px ( top ) sample for the all jets channel. 

The total systematic error is the order of 3.0 GeV for the high pT sample. 
This value is larger than in the case of the lepton plus jets channel where the 
top mass is determined in the same way, as the invariant mass of the three 
jets coming from the hadronic top decay (see section 2). Clearly, the sources of 
systematic uncertainties have an impact on the resolution of the kinematic fit. 

4.5 Summary 

It has been shown that the top mass can be determined in the all jets channel. 
The ti signal is extracted from the huge QCD background {S/B ~ 3 x 10~* at 
production level) by the use of kinematic cuts and a kinematic fit which allows to 
reconstruct the complete final state topology. The signal over background ratio 
can be further increased by selecting high px events. Once the entire tt decay is 
reconstructed, the top mass is determined as the invariant mass of the three jets 
arising from each top quark {t — s- jjb and t — > jjb). It was shown that a total 
error on the top mass of the order of 3 GeV can be reached. 

5 Top mass measurement in leptonic final states with 

Here, one exploits the correlation between the top mass and the invariant mass 
of the system made of a J/'F from the decay of a b-hadron and the isolated 
lepton (e or /i) coming from the associated W decay (see figure ISTjl 32 . In order 
to uniquely define the final state topology and therefore to reduce considerably 
the combinatorial background, the presence of a muon-in-jet (with the same sign 
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Fig. 37. Diagram of the top decay to leptonic final state with J/tf'. 
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than the isolated lepton) from the b-quark is required in the other top quark 
decay. 

The overaU branching ratio is 3.2 x 10~^. Due to this strong suppression this 
method will be only applicable during the high luminosity phase of the LHC, 
where 2700 events will be produced per year. 



5.1 Analysis 

Events are selected by requiring one isolated lepton with pr > 30 GeV and 
\r)\ < 2.4, and three non-isolated muons with px > S GeV and \r]\ < 2.4, with 
the invariant of two of them (with opposite signs) being compatible with the 
J/W mass. After selection cuts are applied and for one year of running at high 
luminosity, about 430 events are expected. The IJ/\P invariant mass distribution, 
for five years high luminosity running, is shown on figure 

The distribution fitted by a Gaussian leads to a mean IJ/'I' invariant mass of 
68.1 GeV with a peak width of 22.4 GeV. The combinatorial background is small. 
For five years of running at high luminosity, the statistical error on M;j/^ would 
be approximately 0.5 GeV. To improve the statistics, various strategies might be 
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considered, such as for example removing the requirement of the muon-in-bjet 
and determine the final state topology by jet charge measurements. Nevertheless, 
further study is needed to evaluate the effectiveness of these approaches. 

The correlation between the top mass and IJ invariant mass is shown on 
figure 1211 One expects the uncertainties of the top mass to scale as a factor 
1/0.58 ~ 1.7 compared to the estimated errors on Mu/^. Therefore, for five 
years of running at high luminosity, the statistical uncertainty on the top mass 
would approximately 0.8-0.9 GeV. 

5.2 Background processes 

The major sources of background come from processes involving hh production. 
Potential backgrounds such as W+jets, Z+jets and boson pair productions are 
briefly discussed. 

The boson pair production processes (WW, WZ and ZZ) have small cross- 
sections compared to signal Furthermore, all the final states topologies are 
different with respect to the signal (except for one particular case but with tiny 
cross-sections |3S1)- These processes can be neglected. 

The cross-sections for W-|-jets and Z-|-jets processes are about a factor of 20 
higher or similar to the signal cross-section, respectively. Here also nevertheless, 
the final states are different with respect to the signal. These processes can be 
neglected as well. 

The total cross-section for inclusive bb production is approximately 10^ larger 
than the signal cross-section. The final state is simply missing one isolated lepton 
compared to the signal. Taking into account the various branching ratios and 
for one year of running at high luminosity, of the order of 3.5 x 10^ && events 
with three non-isolated muons can be expected, compared to 2700 signal events. 
10^66 events with three non-isolated muons have been generated. Only 2.9% of 
the events are reconstructed with three non-isolated muons with pT > i GeV 
and |?7| < 2.4. Only 0.05% of the events are reconstructed with one isolated 
lepton with px > 30 GeV and \ri\ < 2.4. No events survive when both of these 
requirements (corresponding to the signal selection cuts) are applied. In addition, 
a cut on the transverse missing energy could be applied (p™*'*'* > 20 GcV) which 
would reduce the bb rate by a factor of 6 and have no effect on the signal |33| . 
Therefore, and although more statistic would be helpful, the bb background 
process can probably be controlled. 

The Wbb process has exactly the same final topology as the signal. The 
cross-section, for W Iv, is approximately 85 pb (Mj- After selection cuts, the 
reconstruction efficiency if 1.2% (16% for the signal). Therefore, after selection 
cuts, the signal over background ratio would S/B « 55. This background process 
can also be controlled. 

In conclusion, all the background processes are either negligible or can be 
kept under control. 

5.3 Systematic uncertainties 

This technique is insensitive to the jet energy scale which is the main source 
of systematic uncertainty in direct top mass measurements. However, the main 
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Fig. 40. Dependence of the reconstructed Mij/^ on the top quark transverse 
momentum. 



limitation to a precise determination of the top mass using this method relies on 
how well the Monte Carlo describes the top production and decay. Particularly, 
the proper description of the fragmentation of the b hadrons is a crucial point. 

The most relevant sources of systematic uncertainties have been investigated 
in the following. Shifts of M/j/,j/ {AMij/,p) have been determined, defined as 
the difference between the value of Mu/^ under nominal conditions and under 
the condition of the systematic uncertainty source, as described in the previous 
sections. Due to the huge amount of the required Monte Carlo statistics, the 
mass shifts have been obtained with a non negligible error 33 . Nevertheless, the 
quoted numbers should give a realistic estimate of the impact of the systematic 
uncertainty sources. 

The mass shifts obtained for initial state radiation, parton distribution func- 
tion and the b-quark fragmentation parameter are summarized in table IT^ 



Source 


Mass shift in GeV 


ISR 


0.1 


PDF 


0.2 


£& ± 10% 


+0.4/0.0 


tb ± 20% 


+0.8/0.1 



Table 12. Muj^ shifts for the various systematic uncertainties sources. 

The IJ/^ invariant mass can be determined with a systematic uncertainty 
of the order of 0.5 GeV which translates to a systematic error on the top mass 
of the order of 1 GeV. 
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5.4 Top transverse momentum: 

The stability of Mu/^ as a function of pritop) has been controlled and a strong 
Pritop) dependence has been found, as shown on figure HUl This analysis has 
been repeated at generator level using the following selection criteria: no cut 
applied, cuts on isolated lepton only, cuts on non-isolated muons only, and all 
cuts applied. As expected, no dependence is observed when no cuts are applied. 
The cuts on the non-isolated muons introduces a small effect whereas the cuts 
on the isolated lepton has a strong impact PH]- This effect is a kinematic one. It 
can be kept under control as long as the Monte Carlo is well tuned to the data, 
which is needed in any case for the determination of the top mass from Muj^. 

5.5 Summary 

During the high luminosity phase of the LHC, the top quark mass can be deter- 
mined in leptonic final states with J /^. This indirect method relies heavily on 
the proper Monte Carlo description of the top production and decay. The top 
mass can be determined with both a statistical and systematic uncertainty at 
the level of 1 GeV. 

6 Conclusion 

The LHC will be an excellent place to study the top quark properties. The 
very large sample of top events that will be accumulated will allow a precision 
measurement of the top quark mass. Various methods applied to statistically 
independent samples gathering all dominant decay channels of the top quark 
have been investigated. The studies have shown that after only one year of data 
taking at low luminosity (per 10/6~^), a total error on the top mass at the level 
of 2 GeV can be achieved. In the inclusive lepton plus jets channels, the error 
can probably be further reduced down to 1 GeV. In all channels, the errors are 
dominated by systematic uncertainties (the systematic errors are summarized in 
table 

The analyses presented in this paper are differently sensitive to the various 
sources of systematic errors. This will allow reliable cross-checks between the var- 
ious methods and an efficient extraction of the combined ATLAS measurement 
of the top quark mass. 

In the inclusive lepton plus jets channel, errors are dominated by the b-jet 
energy scale and the knowledge of FSR. It was shown that these effects can be 
better controlled using a continuous jet definition. A possibility has been seen to 
reduce the systematic error due to the b-jet scale uncertainty by calibrating the 
b-jets with the same calibration as determined for light jets, though this must be 
further studied will full detector simulation before conclusions can be reached. 
In the lepton plus jets channel, using a sub-sample of hight pt top events, the 
top mass can be reconstructed with a large calorimeter cluster. In the all jets 
channels, it was demonstrated that the tt signal can be efficiently extracted from 
the huge QCD background. In the dilepton channel, it was presented that despite 
the two undetected neutrinos the final state can be fully reconstructed assuming 
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Source of error 


Lepton+jets 


Lepton+jets 


Dilepton 


All jets 


in GeV 


inclusive 


large clusters 




high pT 




sample 


sample 




sample 


Energy scale 










Light jet energy scale 


0.2 


- 


- 


0.8 


b-jet energy scale 


0.7 




0.6 


0.7 


Mass scale calibration 




0.9 






UE estimate 




1.3 






Physics 










Background 


0.1 


0.1 


0.2 


0.4 


b-quark fragmentation 


0.1 


0.3 


0.7 


0.3 


Initial state radiation 


0.1 


0.1 


0.1 


0.4 


Final state radiation 


0.5 


0.1 


0.6 


2.8 


PDF 






1.2 





Table 13. SumniMry of the systemMics errors in the top mass measurement, in 
the lepton plus jets channel, in the all jets channel and in the dilepton channel. 



a value for the top mass. Finally, during the high luminosity phase of the LHC, 
the top mass can also be precisely determined in leptonic final states with J/^. 
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